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1 .  INTRODUCTION 


For  high-power  microwave  emission  from  relativistic  reflex  triodes  and 
diodes  to  be  analyzed  (for  a  uniform  cylindrical  beam  of  radius  R),  we  must 
first  determine  the  average  power  spectral  density  per  steradian  radiated 
perpendicular  to  the  beam.  For  the  relativistic  sheet  model  of  the  electron 
dynamics,  that  density  is  given  by1 

d3Pi/do)d2fl  =  CR'*a)2|s(u,)|2F2(tl)R/c)  (1) 

S(w )  =  I *(w),  from  Brandt  et  al,1  with  C  =  pqtio2/4cT,  and 

S(u>)  =  (1/2r)  /“  dt  exp(iut)  S(t)  .  (2) 


Moreover , 


and 


N(t) 

=  Z  VtJ 

m=1 


(3) 


F2(x)  =  4J^(x)/x2  .  (4) 

In  these  terms,  d2£l  is  the  differential  solid  angle,  yQ  is  the  permeability  of 
free  space,  a  is  the  sheet  surface  charge  density,  R  is  the  beam  radius,  c  is 
the  speed  of  light,  T  is  the  voltage  pulse  duration,  u>  is  the  angular 
frequency,  N(t)  is  the  number  of  sheets  in  the  source  region  at  time  t,  vffl(t) 
is  the  velocity  of  the  m  sheet  at  time  t,  and  is  the  first-order  Bessel 
function  of  the  first  kind.  Equation  (1)  was  obtained  by  solving  the  wave 
equation  with  a  current  distribution  represented  by  the  sheet  dynamics, 
including  only  radiative  time  delays  which  are  perpendicular  to  the  beam, 
while  using  the  dipole  far-field  approximation.  S(t)  is  the  sum  of  the  sheet 
velocities.  This  introduces  the  function  which  is  at  the  heart  of  the  present 
study.  New  let  us  discuss  what  is  to  be  done  with  it  and  why. 


The  purpose  of  this  paper  is  to  investigate  the  extrema  of  equation  ( 1 ) . 
By  extrema  we  mean  maxima,  minima,  or  inflections  in  one  dimension  and,  on 
surfaces,  their  two-dimensional  analogs.  Such  an  understanding  is  necessary 
to  the  optimal  design  of  high -power  microwave  sources. 

It  is  known  that  S(t)  is  a  very  complicated  function  of  time.  S(t)  has 
some  general  features,  however.  It  is  zero  for  t  <  0  and  for  t  >  0  makes  a 
rapid  transition  to  a  regime  in  which  aperiodic  oscillations  about  a  mean 
value  occur  over  a  period  of  time  corresponding  to  many  periods  of 
oscillation.  Ignoring  the  superimposed  aperiodic  fluctuations,  a  simple,  very 
coarse-grained  representation  of  this  behavior  would  be  to  set 

S(t)  -  NavaH(t)  ,  (5) 


1 Howard  E.  Brandt,  Alan  Bromborsky,  Henry  B.  Bruns,  H.  Allen  Kehs,  and 
George  P.  Lasche,  Gigawatt  Microwave  Emission  from  a  Relativistic  Reflex 
Triode,  Harry  Diamond  Laboratories,  HDL-TR-1917  (August  1980). 


p 


where  N&  is  the  mean  number  of  sheets,  va  is  the  mean  velocity,  and  H(t)  is 
the  Heaviside  step  function  given  by  equation  (6)  below.  From  equation  (2) 
one  finds  S(oi)  =  S/w  with  S  =  iNava/2ir,  and  we  are  using  u  as  shorthand  for  u 
+  i6  with  6  >  0+: 


H  (t ) 


0  t  <  0 
1  t  >  0 


(6) 


Inserting  the  Fourier  transform  of  equation  (5)  into  equation  (1)  gives  (K  = 
S2C): 


d3Pi/dud2£J  =  KR^F2  (ojR/c  )  .  (7) 

We  observe  that  in  this  naive  representation  the  u>  dependence  appears  entirely 
in  F2.  One  of  the  principal  problems  just  mentioned  involves  maximizing 
equation  (7)  with  respect  to  u  for  some  given  R.  We  see  that  this  amounts  to 
finding  extrema  for  F2  alone  given  the  very  simple  S(m)  corresponding  to 
equation  (5).  In  the  first  part  of  section  2  we  consider  the  extrema  of  F2(x) 
with  x  =  u)R/c  and  show  that  minina  and  maxima  occur  at  the  zeros  of  J1  and  J2» 
respectively.  Simple  asymptotic  expressions  for  F2,  J1 ,  and  J2  are  given  for 
small  and  large  arguments.  Where  S(to)  cannot  be  considered  to  be  of  the  form 
corresponding  to  equation  (5),  we  must  instead  seek  the  extrema 
of  |  S  (w )  1 2 (<oR/c )  at  constant  R.  This  problem  is  the  subject  of  section 
2.4.  In  section  3,  we  turn  our  attention  to  the  variation  of  equation  (1) 
with  R  at  constant  u.  In  particular,  we  seek  extrema  of  a  function 

f2  =  x**F2(x)  =  4x2J^(x)  ,  (8) 

with  x  *  uSR/c  as  before.  This  function  has  minima  at  J,  *  0  and  maxima  where 
JQ  *  0.  Asymptotic  forms  are  presented  for  f2,  JQ,  and  J.,  for  small  and  large 
arguments.  In  section  4,  we  consider  the  surface  given  by 

h2  -  R^S^^CtoR/c)  ,  (9) 

and  seek  extrema  with  respect  to  u  and  R.  Inspection  of  equation  (9)  reveals 
that  absolute  minima  should  occur  where  S(u>)  or  J1  vanishes.  The  former 
condition  defines  a  set  of  planes  of  constant  ui,  and  the  latter  a  set  of 
hyperbolas  along  which  ojR  is  constant.  The  maxima  are  found  by  solving  two 
simultaneous  algebraic  equations  for  u  and  R.  Clearly,  a  knowledge  of  S(uj)  is 
required  for  a  complete  solution. 

In  the  work  that  follows,  the  available  literature  on  Bessel  functions  has 
been  used  extensively.  The  theory  of  Bessel  functions  is  well  summarized  in 
texts  by  Arfken2  and  Hildebrand.3  Greater  theoretical  detail  is  supplied  by 


2G.  Arfken,  Mathematical  Methods  for  Physicists,  New  York,  Academic  Press 
(1968). 

3F.  Hildebrand,  Advanced  Calculus  for  Applications ,  prentice- Hal 1 : 
Englewood  Cliffs,  NJ  (1962). 
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Watson.4  Frequent  use  has  also  been  made  of  formulae  summaries  to  be  found  in 
mathematical  compendia  edited  by  Gradshteyn  and  Ryzhik5  and  by  Abramowitz  and 
Stegun.6  Also  found  useful  was  the  short  table  of  Bessel  functions,  Jn,  of 

order  0,  1,  and  2  for  values  of  the  arguments  0  <  x  <  17.5  in  steps  of  0.1 

(Stegun).  For  more  accurate  work,  a  reference  of  the  Harvard  University 
Computation  Laboratory  has  been  very  useful:  Tables  of  Bessel  Functions  of 
the  First  Kind,7  for  0  <  x  ^  100  with  incremental  steps  of  0.001  up  to  x  =  25 

and  0.01  thereafter.  Lastly,  Boas8 9  clearly  describes  the  process  of 

determining  extrema  for  unbounded  surfaces  as  well  as  for  those  that  are 
bounded.  The  latter  case  would  occur  if  there  were  lower  or  upper  bounds  of  a 
practical  nature  on  <u  and  R. 


2.  EXTREMA  AT  CONSTANT  BEAM  RADIUS 

2.1  Asymptotic  Forms  of  Interference  Form  Factor  F2  and  Approximate 
Location  of  Minima  and  Maxima 


Figure  1  illustrates  the  behavior  of  F2  for  0  £  x  £  7,  where  F2  is 
given  by  equation  (4).  Its  value  at  x  =  0  is  precisely  1.  The  first  two 
zeros  of  the  function  occur  at  3.83  and  7.02,  J.j(x)  being  zero  here.  In  fact, 
J1  =  0  for  x  *  0  locates  all  the  minima  of  the  function.  Indeed,  one  might 
expect  as  much  by  inspection  of  equation  (4).  The  graph  indicates  two 
successive  maxima.  The  first,  at  x  =  0,  is  much  larger  than  the  second,  which 
occurs  at  x  =  5.14.  This  secondary  peak  has  a  value  of  1.75  x  10~2.  We  will 
show  that  these  peaks  fall  off  as  x“3  for  large  x  and,  moreover,  are  located 
at  the  zeros  of  J2* 

J1  may  be  represented  by  a  few  terms  of  the  series  (n  =  1 ) : 

co  g 

V>  -  W?  £  .Hi"  n).  - 

8=0 

2 

for  small  arguments.  Using  this  we  find  that  F*  is  given  by 
F2  =  1  -  x2/4  +  5x4 / 1 92  -  7x6/4608  +  ...  . 


4G.  S. Watson,  A  Treatise  on  the  Theory  of  Bessel  Functions,  2nd  ed., 
Cambridge,  England,  The  University  Press:  Sew  York,  The  MacMillan  Company 
(1944). 

5 i.  Gradshteyn  and  I.  Ryzhik,  Table  of  Integrals,  Series,  and  Products,  4th 
ed.,  A.  Jeffrey,  ed..  Sew  York,  Academic  Press  (1980). 

6M.  Abramowitz  and  I.  Stegun,  eds..  Handbook  of  Mathematical  Functions  with 
Formulas,  Graphs,  and  Mathematical  Tables,  Washington,  D.  C.,  U.  S.  Government 
Printing  Office  (1964). 

7 Harvard  University  Computation  Laboratory,  Tables  of  the  Bessel  Functions 

of  the  First  Kind  of  Orders,  Cambridge,  Harvard  University  Press  ( 1947-51) . 

9M.  Boas,  Mathematical  Methods  in  the  Physical  Sciences,  Sew  York,  Wiley 
(1966). 


(10) 


(11) 


Structure  Factor  (F2) 


Argument  (x) 


The  first  four  terms  adequately  represent  F  up  to  x  =  1.6  (less  than  1 
percent  error)  as  shown  in  table  1.  The  large  argument  approximation  for  JR 
is 


For  n  *  1 , 

Making  use  of  this. 


Jn  a  (2/wx)1/2  cos  (x  -  nw/2  -  ir/4) 
J1  ■  (2/irx)1/2  sin  (x  -  ir/4)  . 

one  finds  for  F2, 


(12) 

(13) 


F2  =  4(1  -  sin  2x)/irx3 


(14) 


TABLE 

1.  SMALL  ARGUMENT 

OF  FORM  FACTOR 

REPRESENTATION 

F2 

X 

F2  (series) 

F2 (exact) 

Deviation  («) 

0.0 

1. 

1. 

_ 

0.2 

0.9900 

0. 9900 

— 

0.4 

0.9607 

0.9607 

— 

0.6 

0.9133 

0.9133 

— 

o.e 

0. 8503 

0. 8503 

— 

1.0 

0.7745 

0.7746 

0.01 

1.2 

0.6895 

0. 6897 

0.03 

1.4 

0. 5986 

0. 5994 

0.13 

1.6 

0. 5052 

0.5075 

0. 45 

1.8 

0.4117 

0.4175 

1.39 

2.0 

0.3194 

0.3326 

3.97 

This  expression  shows  that  the  successive  maxima  of  F2  for  large  x  fall  off 
rapidly  as  x“3.  Equation  (14)  may  be  used  to  determine  the  large  argument 
locations  of  zeros  (minima)  as  well  as  maxima  of  F2»  The  minima  occur  for  sin 
2x  ■  +1,  requiring 


Xjj  =  x(n  +  1/4) 


(15) 


Zeros  as  predicted  by  equation  (15)  are  compared  with  the  actual  values  as 
determined  by  J,  =  0  and  n  or  x  *  0,  in  table  2.  Results  are  more  than 
adequate  for  n  >  2.  Appendix  A  shows  that  Bessel  functions  of  odd  order,  J1 , 
J3,  etc.,  have  zeros  which  obey  equation  (15)  for  large  x. 

The  maxima  of  F2  occur  where  sin  2x  =  -1 .  Here 


xn  =  x(n  +  3/4)  . 


TABLE  2. 


COMPARISON  OF  APPROXIMATE  FORMULA  WITH 
J1  =  0  IN  PREDICTING  MINIMA  OF  F2  (MINIMA 
ARE  PREDICTED  FOR  N  >  1) 


(16) 


n 

Actual  zero 

Approximate 

zero 

Deviation  (%) 

0 

0 

0.79 

__ 

1 

3.83 

3.93 

2.54 

2 

7.02 

7.07 

0.71 

3 

10. 17 

10.  21 

0.39 

4 

13.32 

13.35 

0.22 

5 

16.47 

16.49 

0.12 

6 

19.62 

19.63 

0.05 

7 

22.76 

22.78 

0.09 

8 

25.90 

25.92 

0.08 

9 

29.05 

29.06 

0.03 

10 

32.19 

32.20 

0.03 

We  compare  the  maxima  locations  as  predicted  by  equation  (16)  with  values 
given  by  setting  J2  =  0,  x  *  0  in  table  3.  There  is  good  agreement 
(percentage  of  deviation  less  than  1  percent)  for  n  >  3.  We  show  in  appendix 
A  that  Bessel  functions  of  even  order,  Jq,  3^,  etc.,  have  zeros  which  obey 
equation  (16)  for  large  x. 


TABLE  3.  ACTUAL  AND  APPROXIMATE  LOCATIONS  OF  F2  MAXIMA 


n 

Actual  location 
<J2  «  0) 

Approximate 

location 

Deviation  (%) 

0 

0 

2.36 

_ 

1 

5.14 

5.50 

7.0 

2 

8.42 

8.64 

2.6 

3 

11.62 

11.78 

1.4 

4 

14.80 

14.92 

0.8 

5 

17.96 

18.06 

0.6 

6 

21.12 

21;  21 

0.4 

7 

24.27 

24.35 

0.3 

8 

27.42 

27.49 

0.3 

9 

30.57 

30.63 

0.2 

10 

33.72 

33.77 

0.1 

9 


Hum 


I 


Inspection  of  equation  (14)  reveals  that  the  nagnitude  of  F2  maxima 
for  large  Xjj  is  given  by 

F2  =  8/irx3  .  (17) 

Equation  (17)  is  compared  with  equation  (4)  in  table  4. 

TABLE  4.  COMPARISON  OF  EXACT  AND  APPROXIMATE  MAGNITUDES  OF 
SUCCESSIVE  F2  MAXIMA 


n 

VJ2  -  °> 

4J^/x2 

8/wx3 

Deviation  ( % ) 

0 

0 

1 

m 

_ 

1 

5.136 

0.017498 

0.018796 

7.42 

2 

8.417 

0.004158 

0.004270 

2.69 

3 

11.620 

0.001601 

0.001623 

1.37 

4 

14.796 

0.000779 

0.000786 

0.90 

5 

17.960 

0.000437 

0.000440 

0.69 

6 

21.117 

0.000269 

0.000270 

0.  37 

7 

24.270 

0.000178 

0.000178 

— 

8 

27.420 

0.000123 

0.000124 

— 

9 

30.570 

0.000089 

0.000089 

— 

2. 2  Exact  Location  of  Maxima  and  Minima  for  F2 

In  determining  maxima  and  minima,  it  will  be  helpful  to  use  the 
recurrence  relations 

Jn-1  =  nJnA  +  JA  (18) 

and 

Jn+1  =  nJn/*  -  JA  *  <19> 

Taking  the  first  derivative  of  F2  and  setting  the  result  to  zero,  we  get 

dF2 /dx  =  aJ^xJJ  -  J^/x3  =  0  . 

From  equation  ( 1 9)  with  n  *  1 ,  we  have 

xJi  -  Ji  =  _xJ2  * 

Thus, 

dF2/dx  -  -Q3^J2/x2  =0  (20) 

is  a  condition  for  relative  extrema. 

The  second  derivative  of  F2  is 


d2F2/dx2  -  ISx-3^^  -  8x"2(J1,J2  +  ^  J£) 


(21) 


Inspection  of  equation  (20)  reveals  that  it  is  satisfied  for  J1  or  J2  =  0, 
provided  x  *  0.  At  x  =  0,  J1  and  J 2,  as  well  as  x,  simultaneously  go  to 
zero.  We  will  investigate  this  situation  last.  Let  us  first  note  that  if 
x  *  0,  J1  and  J2  do  not  have  zeros  in  common.  Setting  J1  =  0,  but  J2  and 
x  *  0  satisfies  equation  (20).  In  equation  (21),  we  have 

d2F2/dx2  =  -8J.{  J2/x2  . 

Setting  n  =*  1  in  equation  (19),  and  noting  that  =  0,  we  have  =  -J2,  so 

d2F2/dx2  =  8x“2J^  >  0  , 

and  there  are  minima  here  (J^  =0).  On  the  other  hand,  setting  J2  -  0, 

d2F2 /dx2  =  -8x~2 J, . 

Now  put  n  =  2  in  equation  (18),  obtaining  J£  =  J1  with  J2  =  0.  Thus,  here 

d2F2/dx2  =  -8x“2jf  <  0  , 

and  we  have  maxima  at  these  x.  Finally,  making  use  of  equation  (11),  we  find 

dF2/dx  =  -x/2  +  0(x3)  , 

which  goes  to  zero  as  x  -►  0,  indicating  a  maximum  or  minimum  here.  Going 
further, 

d2F2/dx2  =  -1/2  +  0(x2 )  , 


so  that  at  x  =  0 


d2F2/dx2  =  -  1/2  <  0  , 

indicating  a  maximum.  At  this  point  F2  =  1,  as  indicated  by  equation  (11). 

Summarizing,  we  have  maxima  at  the  zeros  of  J2  (now  including  x  =  0) 
and  minima  at  the  zeros  of  J1  with  x  *  0. 

2. 3  Frequencies  for  Minima  and  Maxima  of  F2  at  R  «  2.54  cm 


In  general,  the  frequencies,  fn,  at  which  maxima  of  F2  occur  if  R  is 
held  constant  are  given  by 

fn  -  cxn/2wR  ,  (22) 

where  Xj^  is  a  zero  of  J2«  The  largest  value  of  F2  occurs  at  xQ  =  0.  The 
successive  maxima  of  F2  are  given  exactly  by  equation  (17)  and  fall  off  as 
x“3  for  sufficiently  large  Xj,.  For  large  x„  we  have  the  approximate 
relations 

fn  -  c(n  +  3/4)/2R  ,  (23) 
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and 


(24) 


F2  =  8/ir4(n  +  3/4)3 

following  from  equations  (16),  (17),  and  (22). 

The  f  values,  where  F2  =  0,  are  also  given  by  equation  (22),  but  now 
the  ^  are  zeros  of  J,  (x)  for  x  *  0.  For  large  enough  x, 

fn  =  c(n  +  1/4)/2R  .  (25) 

We  summarize  exact  values  of  fn  at  minima  and  maxima  in  table  5  for  R  = 
2.54  cm  (c/2R  =  5.9  GHz).  Exact  values  of  F2  for  the  maxima  will  be  found  in 
table  4.  The  first  two  cycles  of  F2  plotted  versus  frequency  are  presented  in 
figure  2. 


TABLE  5.  FREQUENCIES  (IN  GHz)  FOR  MINIMA  AND  MAXIMA  WITH 
BEAM  RADIUS  =  2.54  cm 


n 

Minina 

Maxina 

Exact 

Approximate 

Exact 

Approximate 

0 

_ 

1.48 

0 

4.43 

1 

7.19 

7.  38 

9.65 

10.33 

2 

13.18 

13.28 

15.81 

16.  23 

3 

19.10 

19.18 

21.82 

22.13 

4 

25.02 

25.08 

27.79 

28.03 

5 

30.93 

30.98 

33.73 

33.93 

6 

36.85 

36.88 

39.66 

39.  83 

7 

42.74 

42.78 

45.58 

45.73 

8 

48.64 

48.68 

51.50 

51.63 

9 

54.56 

54.58 

57.41 

57.53 

F2  =  4J^(x)/x2  versus 


Frequency  (GHz) 
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2.4  Other  Beam  Power  Spectra 


If  S(<d)  *  S/oj, 
extremization  is 


the  appropriate  factor  in  equation  ( 1 )  to  subject  to 


g2  =  |  S (to)  J  J^(uR/c) 


(26) 


By  inspection  of  equation  ( 26)  we  expect  absolute  minima  to  occur  where  either 
|s(u)j  or  J^CuR/c)  =  0.  These  conditions  yield  specific  planes  beginning  with 
u>  =  0  which  are  perpendicular  to  the  u  axis  and  cut  the  surface  prescribed  by 
equation  (1).  Taking  the  first  derivative  of  g2  and  setting  it  equal  to  zero 
gives 

dg2/dcu  =  2|s(a))|j1[J1|s((l))|'  +  R  j  S  (a»)  |  J^/c]  =  0  ,  (27) 

where  |s(u>)| 1  and  Jj  denote  derivatives  of  |s(w)|  or  J1  with  respect  to  their 
arguments.  The  factor  in  parentheses  yields  extrema  of  equation  (26)  other 
them  absolute  minima  (i.e.,  extrema  where  g2  is  nonzero).  Among  these  may  be 
maxima  of  varying  height  which  include  one  or  more  absolute  maxima,  other 
minima  (located  between  the  peaks)  which  give  nonzero  g2,  as  well  as  points  of 
inflection.  These  cases  are  sorted  out  by  taking  the  second  derivative. 
Making  use  of  equation  (19)  we  have 


J|  —  cJ,  /o>R  —  J2  . 

On  substituting  the  above  expression  in  equation  (27),  we  get  as  the  condition 
for  extrema  (barring  absolute  minima) 

[|s(w)|'  +  oi-1|s(o))|]ji  -  r|s(«)|j2/c  =  0  .  (28) 

Note  that  equation  (28)  reduces  to  the  earlier  condition  for  maxima,  namely, 
J2  *  0,  if  |S(w)|"  =  S/w  The  test  for  maxima  is  given  in  appendix  B,  where 
it  is  shown  that  the  behavior  of  | S ( u> )  |  " /S ( oj )  is  the  determining  factor. 


3.  EXTREMA  AT  CONSTANT  FREQUENCY 

3.1  Approximate  Representations  of  Form  Factor  f2 


Figure  3  illustrates  the  behavior  of  f2  (eq  (8))  for  0  <  x  <  7.  It 
has  two  peaks  in  the  interval.  These  are  located  at  2.41  and  5.52,  JQ  being 
zero  at  these  points.  Jq  =  0  also  locates  all  the  other  naxima  of  this 
function.  Also,  as  one  might  expect,  J1  =  0  locates  all  the  minima.  Three  of 
these  zeros  are  seen  on  the  graph.  They  lie  at  x  values  of  0,  3.83,  and  7.02. 

The  maxima  are  seen  to  grow  with  x.  The  first  maximum  has  a 
magnitude  of  6.23;  the  second,  14.11.  Successive  maxima  which  are  off  the 
graph  are  larger  still.  We  will  shew  shortly  that  this  growth  is  proportional 
to  x. 
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Argummt  |«| 


Structure  factor  f2  =  4x2J^(x)  versus 
*  UjR/c. 


From  equations  (8)  and  (11),  we  immediately  determine  a  useful  series 
which  represents  f2  for  small  arguments: 


f2  =  x4  -  x6/4  +  5x8/192  -  7x10/4608  +  ...  .  (29) 


One  obtains  reasonable  values  (less  than  1  percent  error)  up  to  x  =  1.6  using 
only  the  first  four  terms  given  above.  This  is  illustrated  in  table  6. 

TABLE  6.  SMALL  ARGUMENT  REPRESENTATION  OF  FORM 
FACT®  f2 


X 

(aeries) 

f2(«*aet) 

Deviation  («) 

0 

0 

0 

— 

0.2 

0.0016 

0.0016 

— 

0.4 

0.0246 

0.0246 

— 

0.6 

0.1184 

0.1184 

— 

0.8 

0. 3483 

0. 3483 

— 

1.0 

0.7745 

0. 7746 

0.01 

1.2 

1.4291 

1.4302 

0.03 

1.4 

2.2496 

2.3027 

0.13 

1.6 

3.  3108 

3.3258 

0.45 

1.8 

4.3219 

4. 3826 

1.39 

2.0 

5.1111 

5.3218 

3.96 

Substituting  the  large  argument  approximation  for  J1  (eq  (13)  in  eq 
(8))  one  finds 


f2  -  (4x/w)(1  -  sin  2x )  .  (30) 

The  location  of  minima  and  maxima  are  again  given  by  equations  (15)  and  (16), 
respectively.  Table  2  gives  locations  of  the  absolute  minima,  including  n  or 
x  *  0  where  J1  is  zero.  Table  7  gives  maxima  locations  for  f2.  Agreement 
with  the  approximate  formula  and  exact  locations  as  given  by  JQ(x)  -  0  is  very 
good. 
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TABLE  7.  ACTUAL  AND  APPROXIMATE  LOCATIONS  OF  f2  MAXIMA 


n 

Actual  location 
<J0  -  0) 

Approximate 

location 

Deviation  (%) 

0 

2.41 

2.  36 

2.07 

1 

5.52 

5.50 

0.36 

2 

8.65 

8.64 

0.12 

3 

11.79 

11.78 

0.08 

4 

14.93 

14.92 

0.  07 

5 

18.07 

18.06 

0.06 

6 

21.21 

21.21 

— 

7 

24.35 

24.35 

— 

8 

27.49 

27.49 

— 

9 

30.63 

30.63 

— 

10 

33.78 

33.77 

Inspection  of  equation  (30)  shows  that  the  growth  of  maxima  for  f2  is 
given  approximately  by  the  linear  relation 

€n  “  ^n/*  *  (31) 

Equation  (31)  is  compared  with  equation  (8)  in  table  8.  Table  8  also  gives 
local  R  values  for  maxima.  There  are  practical  limitations  on  the  largeness 
of  R.  Figure  4  gives  the  first  three  cycles  of  f2  versus  R  for  35  GHz. 


TABLE  8.  COMPARISON  OF  EXACT  AND  APPROXIMATE  MAGNITUDES 
OF  SUCCESSIVE  f2  MAXIMA 


n 

x(J„  -  0) 

R  (cm) 

4x2J^ 

8x/» 

Deviation  («) 

0 

2.405 

0.33 

6.23 

6.12 

1.77 

1 

5.520 

0.75 

14.11 

14.06 

0.35 

2 

8.654 

1.18 

22.07 

22.04 

0.14 

3 

11.792 

1.61 

30. 05 

30.03 

0.  07 

4 

14.931 

2.04 

38.04 

38.02 

0.05 

5 

18.071 

2.47 

46.04 

46.02 

— 

6 

21.212 

2.89 

54.03 

54.02 

— 

7 

24. 352 

3.32 

62.03 

62.01 

— 

8 

27.49 

3.75 

70.02 

70.00 

— 

9 

30.63 

4.18 

78. 02 

78.00 

— 

3. 2  Exact  Location  of  Maxima  and  Minima  for  Form  Factor  f2 


This  section  gives  the  analytical  method  for  determining  exact  maxima 
and  minima  for  this  case.  The  condition  for  extrema  is 

df2/dx  -  8xJ^( +  xJ})  -  0  .  (32) 
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Structure  Factor  |P) 


From  equation  (18) 


Therefore, 


xJ0  =  J1  +  xJj 


df2/dx  =  SxZJqJ,  =  0  . 

Thus,  extrema  occur  where  JQ  or  J1  is  zero.  Continuing, 

d2f2/dx2  *  16XJQJ,  +  8x2(j^)J1  +  JQJ')  . 


When  J1  *  0, 


d2f 2/dx2  =  8x2JqJ j  . 

From  equation  ( 1 8)  with  J1  =  0,  we  have  J}  =  JQ,  so 

d2f2/dx2  =  8x2J^  >  0, 

and  there  are  minima  here.  On  the  other  hand,  when  JQ  =  0, 

d2f2/dx2  -  8x2J^J1  , 


Figure  4.  Structure  factor 
f2  *  4x2J^(x)  versus 
beam  radius  R  =  cx/u>* 


(33) 

(34) 

(35) 

(36) 

(37) 
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(38) 


but  *  -J1  by  equation  (19)  and 

d^f2/dx2  =  -8x2 <  0  • 

This  designates  the  location  of  maxima. 


4.  SURFACE  EXTREMA 


To  determine  surface  extrema,  we  start  with  equation  (9). 
derivatives  with  respect  to  R  yields 


Taking  partial 


3h2/3R  =  2R|s{ai)  +  uRjj/c)  =  0 


(39) 


Here  denotes  differentiation  with  respect  to  the  argument.  Making  use  of 
equation  (18),  we  have 


3h2/3R  =  2u)R2j1  |s(o))|2j0/c  *  0  .  (40) 

On  the  other  hand,  the  partial  derivative  with  respect  to  u  gives 

3h2/3o)  =  2R2|s(a)))j1(J1d|s(M)|/dto  +  r|s(u>)  |  j*/c)  =  0  .  (41) 

With  the  help  of  equation  (19)  this  is  rendered  in  the  form 

3h2/3<o  =  2R2|s(a))|j1{(djs(a))|/do)  +  |s(aJ)|/w)J1  -  r|s(w)|j2/c}  =  0  .  (42) 

The  bracketed  term  is  identical  to  equation  (28).  As  expected,  absolute 
minima  occur  where  |s(ui)|  or  J1  *  0.  The  former  condition  gives  planes  of 
constant  to,  and  the  latter  yields  a  family  of  hyperbolas  of  constant  ojR. 
Minima  also  occur  where  u  or  R  is  zero.  Inspection  of  equations  (40)  and  (42) 
reveals  that  other  extrema  (not  including  the  absolute  minima)  occur  if  we 
satisfy 


J0  * 


and 


(43) 


(d|s(u))|/du>  +  |  S  (u>)  |  /o»)  J,  -  r|s(u))|j2/c  *  0  (44) 

simultaneously.  Let  the  zeros  of  be  denoted  new  by  xnQ.  They  are  given  in 
table  8.  Then 


R  *  cxn0/u>  (45) 

replaces  R  in  equation  (44).  Do  not  forget  the  arguments  of  J1  and  J2I  Thus, 

J-jn  *  Ji(*no)  an<*  J2n  “  J2(*no)  are  now  indeP®ndent  of  u>  and  R.  Equation  (44) 
reduces  to 
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(d|s(u)  |di»)  +  |s(<d)|/<d)J1n  -  xn0js(u)|  J^/u  -  0  .  (46) 

The  solution  depends  on  the  fora  of  |s(o>)|*  Consider  our  earlier  example.  If 
|  S  ((»)  |  »  S/o>»  equation  (46)  reduces  to 


This  condition  can  only  be  met  approximately  for  large  x  where  the  zeroe  of  Jq 
and  J2  converge.  Here  the  extrema,  including  any  absolute  maxima,  occur  along 
hyperbolas  of  constant  wR. 


5.  CONCLUSIONS 

The  problem  of  extremizing  the  average  power  spectral  density  per 
steradian  may  be  approached  from  three  points  of  view.  If  the  operating 
frequency  is  fixed  by  practical  considerations  (i.e.,  the  atmospheric  window 
for  microwave  emission  at  35  GHz),  we  seek  to  vary  R  at  constant  u.  If  R  and 
voltage  are  fixed  by  certain  design  considerations,  but  to  is  not,  we  seek 
extrema  with  varying  co.  When  both  variables  are  unrestricted,  we  must  seek 
extrema  over  the  entire  surface. 

The  (d  dependence  of  the  Fourier  transform  of  the  sum  over  sheet  velocities 
plays  an  important  role  in  the  determination  of  extrema  if  o>  is  not  fixed. 
The  only  exceptions  to  this  are  the  absolute  minima  given  by  the  zeros  of 
J1 .  The  average  power  radiated  is  zero  here.  If  |  S  (&> )  |  is  of  the  form  S/u> 
with  S  some  positive  constant,  maxima  are  located  at  the  zeros  of  J2-  The 
maxima  are  then  seen  to  fall  off  as  u~3  for  large  u>.  This  is  in  addition  to 

the  absolute  minima  associated  with  J1  as  discussed  above.  For  general 

|s(u))|,  additional  absolute  minima  occur  where  |s(<d)|  =0,  the  average  power 
radiated  being  zero  here. 

If  u  is  fixed  and  R  varied,  one  finds  the  absolute  minima  to  be  again 
associated  with  J1  as  expected.  These  minima  are  now  located  at  the 
intersection  of  the  surface,  the  plane  of  constant  u>,  and  planes  of  constant 
R.  Maxima,  however,  are  now  associated  with  the  zeroe  of  JQ.  This  yields 

another  set  of  planes  of  constant  R  associated  with  the  maxima.  The  growth  of 

these  maxima  is  proportional  to  R  for  large  R. 


Allowing  both  u>  and  R  to  vary  freely  leads  to  a  set  of  two  algebraic 


equations  which  must  be  solved  simultaneously  for  u  and  R.  Again,  |S((d) 
J1  ■  0  are  associated  with  absolute  minima.  The  zeros,  if  any,  of  |S(h>) 
planes  of  constant  o>.  On  the  other  hand,  J 


and 

are 


1 


0  now  implies  a  set  of 
hyperbolas  (actually  sheets  bent  in  this  form)  with  u>R  constant.  Their 
intersection  with  the  surface  locates  these  absolute  minima.  As  the  argument 
of  J.,  is  (dR/c,  the  planes  u  or  R  •  0  are  also  associated  with  absolute 
minima.  other  extrema  including  absolute  maxima  are  determined  by  solving 
simultaneously  taro  algebraic  equations  for  the  points  (u,  R)  at  which  the 
extrema  are  located.  For  the  special  case  |  S (<d )  |  «  u)_1S,  the  points  trend  to 
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a  set  of  hyperbolas  which  are  the  simultaneous  solutions  of  JQ  =  J2  =  0  (valid 
only  for  large  arguments).  For  other  |s(ui)|,  the  results  can  be  more 
complicated.  The  importance  of  S(w)  in  determining  the  location  of  maxima 
suggests  the  need  to  develop  more  realistic  expressions  for  theoretical  use  of 
computer  data. 
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APPENDIX  A. — LARGE  ARGUMENT  ZEROS  OF  BESSEL  FUNCTIONS 


Equation  (12)  in  the  body  of  the  report  gives  the  general  expression  for 
Jn(x)  for  large  x.  For  n  =  2s,  we  have  on  substitution 

J2S  =  (2/xx)1/2  cos  ( (x  -  s/4]  -  sir) 

=  (-1 )s(2/xx)1/2  cos (x  -  s/4) 

The  argument  of  the  cosine  factor  is  the  same  for  all  the  J2s*  It  alone 
determines  the  zeros  which  occur  for 

x  =  s(n  +  3/4)  .  (A-1 ) 

On  the  other  hand,  for  odd  n  =  2s  +  1 , 

J2s+1  =  (2/irx)1/2  cos((x  -  s/4]  -  1/2(2s  +  1  )s ) 

*  (-1)s(2/*x)l/2  sin(x  -  s/4)  . 

Therefore,  the  zeros  of  the  J2s+i  whic*1  depend  on  the  argument  of  the  sine 
alone  occur  at 

x  =  s(n  +  1/4)  .  (A-2) 

The  foregoing  remarks  apply  to  sufficiently  large  arguments  of  the  Jn» 
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APPENDIX  B. — CONDITION  FOR  MAXIMA  OR  MINIMA 


As  J.,  and  |s(d))|  in  equation  (19)  *  0,  except  for  absolute  minima,  we  may 
write  When  away  from  these  points 

dg2/d<o  =  2g2  |dln|  S(io)  |  /d<o  +  u"1  -  RJj/cJ^  =  0  ,  (B-1  ) 

and  we  explore  the  region  where  the  last  factor  on  the  right  vanishes.  Talcing 
the  second  derivative  we  have 


d2g2/du2  =  2g2(dln|s(u>)|/du>  +  w"1  -  RJ  /c^)'  .  (B-2) 

The  prime  on  the  last  factor  denotes  differentiation  with  respect  to  to.  As 
2g2  >  0,  whether  we  have  minima  or  naxima  here  depends  on  the  last  factor  of 
equation  (2)  in  the  body  of  the  report.  Talcing  the  indicated  derivative,  we 
have 

g~2d2g2/2du)2  =  |s(u))|"/|s(m)J  -  (  |s(<o)|  v|s(co)|  )2 

-  (o~2  -  (R/c  )2(  J£Jj  -  J^/J,2  .  (B-3) 

We  now  make  use  of  certain  relations  to  simplify  equation  (B-3).  From  the 
vanishing  of  equation  (B-1), 

|s(u>)j '/|s(w)|  =  -w'1  +  RJ 2/c J i  •  (B-4) 

Also,  we  have  the  recurrence  relations  from  equations  (18)  and  (19): 


and 


-x’1  J2  +  J1 


x~lj1  -  J2 


(B-5) 

(B-6) 


On  substitution  of  equations  (B-4)  through  (B-6),  equation  (B-3)  becomes 

g“2d2g2 /2dw2  *  js((i>)|  "/|s(«)|  -  2u-2(l  -  x^/J,)2  -  (R/c)2  .  (B-7) 

This  is  as  far  as  we  can  go  without  a  specific  form  of  |s(<o)|.  Note  that  the 
last  two  terms  are  negative.  Extrema  are  therefore  guaranteed  if 

|s(u)|  "/js(<o)|  <  2to“2(  1  -  xJj/J,)2  +  (R/c)2  .  (B-8) 
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GLOSSARY  OF  TERMS 


c  speed  of  light 

d  distance  between  grid  and  ground  plate 
f  quantity  defined  by  equation  (8) 
g  quantity  defined  by  equation  (26) 
h  quantity  defined  by  equation  (9) 

m  summation  index  introduced  in  equation  (3) 
n  denotes  order  of  Bessel  function 
s  summation  index  in  equation  (10) 

v  sheet  velocity  (v  )  or  average  sheet  velocity  (v_ ) 

m  ci 

x  o)R/c 

C  constant  in  equation  ( 1 ) 

F  form  factor 
H  Heaviside  step  function 

J  Bessel  function  of  first  kind  (nth  order  denoted  by  JR) 

K  constant  in  equation  (7) 

N  upper  limit  on  sum  (equation  (3))  expressed  as  a  function  of  time  or  mean 
number  of  sheets  when  subscripted  by  a 
O  designates  "order  of" 

P  power  or  point  (us,  R) 

R  beam  radius 

5  sum  of  sheet  velocities,  S(t),  or  a  constant 

6  infinitesimally  small  quantity 
1  permeability  of  free  space 

us  angular  frequency 
ft  solid  angle 
x  pi 

o  sheet  surface  charge  density 
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